Cooling many particles to very low temperatures 
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In a recent paper [Beige, Knight, and Vitiello, |quant-ph / 0404160] , we showed that a large number 
of particles can be cooled very efficiently using a bichromatic interaction. The particles should 
be excited by red-detuned laser fields while coupling to the quantized field mode inside a resonant 
and leaky optical cavity. When the coupling constants are for all particles the same, a collective 
behavior can be generated and the cooling rate can be as large as \/N times the single-particle 
coupling constants. Here we study the algebraic structure of the dynamics and the origin of the 
collective cooling process in more detail. 

PACS numbers: 03.67.-a, 42.50.Lc 



I. INTRODUCTION 

Collective features in the microscopic dynamics often 
lead to the emergence of surprising and unexpected ef- 
fects in the evolution of a physical system at the macro- 
scopic level Q. The collective cooling of many two- level 
particles to very low temperatures is discussed here as an 
example of such a macroscopic manifestation of micro- 
scopic collective behavior 2] . It is shown that the collec- 
tive behavior of a large number of particles can produce 
much higher cooling rates than they could be achieved 
by means of individual cooling based on the spontaneous 
decay of the individual particles 3, 4, 5]. 

As in laser sideband cooling techniques for single two- 
level atoms [3, we consider an experimental setup, where 
red-detuned laser fields increase the excitation of the 
particles, thereby continuously reducing the number of 
phonons. Afterwards, the phonon energy can be removed 
constantly from the system. This requires energy dissi- 
pation and yields an overall decrease of the von Neumann 
entropy in the setup 0. One possible decay channel is 
spontaneous emission from the excited states of the par- 
ticles. During such a photon emission, a particle returns 
most likely into its ground state without regaining the 
phonon energy lost in the excitation process. The net 
result is a conversion of the phonons, originally existing 
in the setup in the form of thermal energy, into photons 
escaping the system. 

Here we are interested in realizing much higher cooling 
rates than could be achieved with the help of the above 
described spontaneous decay of individual particles. This 
is possible, when the time evolution of the system remains 
restricted onto a highly symmetric and strongly coupling 
subspace of states throughout the whole cooling process. 
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FIG. 1: (a) Experimental setup for the collective cooling of 
many particles. The particles should be trapped inside an 
optical cavity, where they arrange themselves in the antin- 
odes of the resonator field, (b) Level configuration of a single 
two-level particle for the cooling of a vibrational mode with 
frequency v driven by a red-detuned laser field with Rabi 
frequency 57^ and coupling to a resonant optical cavity with 
strength g. 



Maximum cooling is obtained when the particles exhibit 
cooperative behavior in the excitation step as well as in 
the de-excitation step. To achieve this we assume as in 
Ref. 2] , that the outward energy dissipation is conducted 
by an optical cavity. As shown in Fig. Q the particles 
should be placed inside a resonant optical cavity with a 
relatively large rate for the leakage of photons through 
the cavity mirrors. The particles can then transfer their 
excitation collectively into the resonator field, from where 
the energy leaves the system without affecting the num- 
ber of phonons in the setup and without changing the 
symmetry of the involved states. 

Collective behavior of the system requires furthermore 
that the Rabi frequency Vly of the laser field for the cool- 
ing of a vibrational mode with frequency v and the cav- 
ity coupling constant g are for all particles (practically) 
the same. Initially, the particles should all be prepared 
in their ground state (in the large N limit fluctuations 
can be neglected). As shown below, the collective states 
of the assembly of N particles then experience a very 
strong coupling to the laser field as well as to the cavity 
mode 0. As a consequence, the number of phonons de- 
creases exponentially with a with as large as \/N times 
the single-particle coupling constants, which can be much 
higher than previously predicted cooling rates in compa- 
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rable setups 0, i, H IH [H El El El 13 . 

The proposed coohng scheme might be used to cool a 
large number of particles very efficiently. So it should 
be applicable to the preparation of Bose Einstein con- 
densates. Currently, these experiments mainly use evap- 
orative cooling [13 which systematically removes those 
atoms with a relatively high temperature from the trap. 
Consequently, only a small percentage of the initially 
trapped atoms is finally included in the condensate. If 
one could instead cool all the atoms efficiently, yet at 
the same time avoid the loss of particles, it should be- 
come easier to experiment with large condensates. Cool- 
ing is also crucial for ion trap quantum computing, where 
the achievable gate operation times can depend primar- 
ily on the efficiency of the cooling of a common vibra- 
tional mode E3- First cavity-cooling experiments in- 
volving many particles and observing enhanced cooling 
rates have already been performed 0, 0| . 

In this paper, we consider some algebraic features of 
the dynamics ruling the collective cooling of many par- 
ticles and further clarify some aspects of the underlying 
mechanisms, which have already been studied in Ref. 
A remarkable feature is, for example, the presence of a 
relatively large and negative coherence k^. Instead of 
solving the time evolution of the system explicitly, we 
avoid certain approximations made in Ref. ^ by referring 
to the Heisenberg picture. The system considered in this 
paper and its collective dynamics may be a paradigmatic 
example for other applications of physical interest, like 
the evolution of a system undergoing a continuous phase 
transition, yet preserving some specific features during 
such an evolution. 



II. THE EXPERIMENTAL SETUP 

The cavity coupling constant g is for all particles 
the same, when the particles distribute themselves in 
the antinodes of the resonator field, as discussed in 
Ref. 0, 0| (see Fig. H^a)). Alternatively, a ring res- 
onator could be used, as proposed in Refs. [T3. ll^L ll^ . 
In the following we consider a collection of N two-level 
particles (atoms, ions or molecules) with ground states 
|0)i and excited states (see Fig. Qlb)). The setup 
should be operated in a parameter regime, where 

1/2 

K^VNg, i\/iV?7f7>r with n=(^^nl^ . (1) 

Here k is the decay rate of a single photon in the cavity 
mode, r is the spontaneous decay rate of a particle in 
the excited state and 77 denotes the Lamb-Dicke param- 
eter characterizing the steepness of the trap |2l|. As in 
Ref. 0, we discuss in this paper the two extreme cases, 
namely the cooling of common vibrational modes and the 
cooling of the individual phonon modes in the absence of 
common vibrational modes. 



A. Cooling of common modes 

In the following, denotes the annihilation operator 
for a phonon in the common vibrational mode with fre- 
quency V, while c is the annihilation operator for the 
cavity photons and Ci = |0)ii(l| is the lowering operator 
for particle i. Using the rotating wave approximation |2ll | 
and going over to the interaction picture with respect to 
the interaction-free Hamiltonian, the time evolution of 
the system can be described by the Hamiltonian 

One way to simplify this Hamiltonian is to introduce the 
effective phonon annihilation operator 

h = Y^{^^l^)K with [6,6t] = i^ (3) 
which allows to write the Hamiltonian ^ as 

^^comm = (i77^^CT|6 + 5CTjc-f H.c.) . (4) 



B. Cooling of individual particles 

We now consider the operator hi,^i for the annihilation 
of the phonons of mode v of particle i . Again we assume 
that the corresponding laser Rabi frequencies are for 
all particles the same. Then the Hamiltonian of the sys- 
tem equals in the interaction picture and in the rotating 
wave approximation 

i?indi = (5'?^"^ '^'^i^^-^ + ff^lc + H.c.) . (5) 

Proceeding as above, introducing the effective phonon 
annihilation operator 

h, = ^[^,l^)K,, with {h,h\\ = \ (6) 

and using Eq. the Hamiltonian Q becomes 

ifind = ft^(i?7r!a|6,+5(T|c-}-H.c.) . (7) 

■i 

Although this Hamiltonian has some similarities with the 
Hamiltonian Q , it describes a physically different situa- 
tion. Instead of coupling to a set of common vibrational 
modes, each particle sees his own set of phonons. 

C. Spontaneous emission 

Spontaneous emission is described in the following by 
the master equation |22| 

p = — - p] + K [cpc" — \(^cp — \pc^c) 

+r ^ {(y^p(J\ - \o\aip - ^pajai) (8) 
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with Hi being i?comm or ffi„d, respectively, k being the 
decay rate of a photon in the cavity and F being the de- 
cay rate of the particle excited state. From this equation 
one can easily see that the dissipation of cavity photons 
reduces the energy in the system without affecting the 
state of the particles. Cavity decay therefore does not 
disturb the collective behavior of the system. However, 
this does not apply to the emission of photons from the 
particles, which can negatively interfere with the collec- 
tive cooling process. 



III. BOSONIC SUBSYSTEMS 

In case of the cooling of common vibrational modes, 
the setup consists effectively of two different subsystems. 
One is the two-level particle system. As an effect of the 
emergence of the cooperative behavior, the collection of 
particles manifests itself as a bosonic system. We will see 
below that a further consequence of this is the transition 
to a strong coupling regime, which in turn implies a much 
shorter time scale for the system evolution. The other 
subsystem is the bosonic system of phonons and photons 
with a continuous conversion of phonons into photons. 
We will see that it is convenient to consider a boson mode 
which is a superposition of them. This is analogous to 
the field-atom polariton of Hopfield [2^ 01 ■ 



A. Bosonic behavior of the particles 

Suppose the time evolution of the particles is solely 
governed by the operators (t+ = J2i'^t^ — 
and (73 = - \^)u{0\)- The a, are the Pauh 

matrices. These operators obey the SU (2) commutation 
relations of a fermion-\ike A^-body system 

[0-3,0-=^] = ±CT^ and [ct", cr+] = -2ct3 . (9) 

Under this condition, the time evolution of the system 
with all particles initially prepared in |0)i remains re- 
stricted to a subspace of highly symmetric particle Dicke 
states (for more details see |2,£.5j 26, 27]). Since we as- 
sume a large number A^ of particles, small fluctuations in 
the initial state of the system can be neglected. Let us 
denote by I the number of particles excited by the laser 
action into the upper state |1). Then it can be shown 
that ioi N ^ I the particle system can be described by 
the collective operators 

^^^^T.^^' ^"^;^E^r, ^3^-3, (10) 

with (73 = S+S- - ^N. In the large N limit {N ■» I), 
the operators in Eq. (|10|l obey the relations 

[S3,S^]=±S^, [5'-,5'+] = l. (11) 



This means, the su(2) algebra © written in terms of 
and 5*3 contracts in the large- A^ limit to the (projective) 
e(2) (or Heisenberg-Weyl) algebra ijTTjl pll29ll30l|. 

The physical implication of the familiar commutator 
relations (|11|1 is that the particles behave no longer like 
individual fermions. Instead, the time evolution gener- 
ates the excitation of bosonic modes, namely collective 
dipole waves, with denoting the creation and anni- 
hilation operators of the associated quanta obeying the 
usual commutation relation (|ll|l . As a consequence, the 
collection of the single two-level particles manifests it- 
self as a bosonic system. In other words, the ladder of 
equally-spaced Dicke states approximates to a weakly- 
excited harmonic oscillator. 

Using the operators (|10|) . the Hamiltonian Q for the 
cooling of common vibrational modes can simply be writ- 
ten as 

iJcomm = h{x S+b + y S+C + B..C.) , (12) 

where we have introduced the notation 

x=^\^rjn and y = y/N g . (13) 

From this one sees that the time evolution of the system 
is mainly governed by the parameters x, y and k, which 
scale as VA- We thus have, as a consequence of the 
emergence of the particle collective behavior, the tran- 
sition to the strong coupling regime, ^rj^ ^VNrjfl 

and g —> \/Ng. The evolution of the system happens no 
longer on the time scale given by the parameters ■^rjil and 
g but on a much shorter time scale defined by ^y/Nrjfl 
and VNg. 

B. Photon-phonon exchange 

Another factor that contributes significantly to the de- 
scribed cooling mechanism is the interface or ambivalent 
role played by the particles with respect to the phonons 
and the photons in the setup. The only difference be- 
tween the particle-phonon and particle-photon interac- 
tion is the difference of the coupling constants, given by 
\/N-^7]fl and \/Ng, respectively. In some sense, the par- 
ticles act as an engine transforming phonon energy into 
photon energy. This last one is then dissipated outward 
through the coupling with the cavity which allows en- 
ergy leakage. However, in the absence of the leakage of 
photons through the cavity, the inverse transformation, 
photons to phonons, is also possible in principle. In such 
a situation, interesting interference and coherence effects 
arise, which we analyze below in detail. 

Indeed, a closer look at Eq. H12() reveals that the Hamil- 
tonian for the cooling of common vibrational modes can 
alternatively be written as 

-ffcomm = hzS+a + R.c. (14) 

with 

z = \/ x'^ + y'^ , a=—{xb + yc) and [a,a^] = l. (15) 
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Instead of interacting with the phonons and photons sep- 
arately, the particles see the boson mode with annihila- 
tion operator a and number operator 



a' a — 



— [x^ 6^6 -t- c'c + yx 



with 



k-i = b''c + bc'^ . 



(16) 



(17) 



Physically, the creation of bosons corresponding to 
does not only affect the number of phonons and the num- 
ber of photons in the system. Inevitably, it also cre- 
ates a coherence between the b and the c subsystem. 
This coherence k^ provides a "symmetric" channel for 
the phonon-photon energy transformation. However, the 
leakage of energy outside the cavity perturbs such a sym- 
metric phonon-photon balancing due to the k^ action. 
The system reacts by subsequent adjustments, trying to 
recover its lost balance. Crucial to such a re-adjustment 
mechanism is the difference in the time scales. The time 
scale for reaching the quasi-stationary state is of the or- 
der 1/iV, while the time scale for spontaneous decay is, 
as we see in the next Section, of order 1 / ^/N. 



IV. COLLECTIVE COOLING OF COMMON 
VIBRATIONAL MODES 

The time evolution of the system turns out to be highly 
non-linear since the second and higher order derivatives 
of the physical observables are much larger than their first 
order derivatives. The system reaches a quasi-stationary 
state on a time scale of the order 1/N. The formation of 
this local equilibrium is solely governed by the effect of 
the Hamiltonian iJcomm- Let us therefore first consider 
the situation, where we can neglect spontaneous emission 
and /v « and F « 0. 



A. Conservation laws in the absence of dissipation 

In this case, the time evolution of the system, governed 
by the Hamiltonian (|14|l . results in a redistribution of 
population between the bosonic subsystem described by 
a^a and the particles described by S^S~. To analyse 
this process we introduce the operators 

Li EE ^{S+a + S-a^), 

Ls = i(5+S- -a^a) (18) 
with the familiar SU (2) commutators 

[Li,Lj] = ie^fc Lfc . (19) 
With this notation, the Hamiltonian (|14|l becomes 



and this, formally, simply generates a rotation around the 
1-axis. In such an algebraic picture, we can immediately 
conclude that the angular momentum L\ and the total 
angular momentum 



I? = L\ 



Li 



= \{S+S- +a'a)[l + \{S+S- +a'a)\ (21) 

are conserved during the time evolution of the system 
under the considered conditions. Especially, the conser- 
vation of implies the conservation of the total num- 
ber of bosons, as accounted for by the number operator 
S+S- +a''a. 

When the number of particles in the excited state |l)i 
remains small compared to TV, the commutator relations 
(|ll|l hold. In such a case, we can assume that the ex- 
pectation value of S^S~ = ^t'^T remains small. 
Neglecting in Eq. (|2U imphes 



(22) 



Using the definition (|15|) , this conservation law translates 
into 



f^{y''b''b + x'c^c~xyk:i) =0, 



(23) 



which coincides with Eq. (16) in Ref. |3| and describes 
the continual balance between the rate of change of the 
expectation value of k^ and the rate of change of the total 
number of phonons and photons in the system. 

To show that Eq. H23|l does not violate the conservation 
of the number of particles in the setup, we remark that 
the system obeys a second conservation law. Considering 
again the Heisenberg picture, we indeed find 



A(6t6 + ctc) = -^(5+5-) = -2zL2 



(24) 



In the presence of many particles, with most of them 
remaining in their ground state, the expectation value of 
S^S^ remains negligible. This implies ■^{S^S~) — 
and, consequently, also L2 = IsJ- Inserting this into 
Eq. H24() , we obtain the particle number conservation law 



= 0. 



(25) 



Hr, 



2hzLi 



(20) 



Despite being a coherence, the expectation value of k^ 
acts like a population. This allows the system to obey 
conservation of the particle number a' a as well as con- 
servation of -|- c^c by balancing the coherence fcs ac- 
cordingly. In the following, we study the conversion of 
phonons into cavity photons and creation of a non-zero 
coherence k^, in more detail. 



B. Quasi-stationary states 

In the absence of spontaneous emission, the system 
reaches after a very short time a stationary state with 
constant expectation values for 6^6, c^c and ^3. To calcu- 
late the corresponding values of the cavity photon num- 
ber and coherence k^ as a function of the phonon number. 
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we use again the Heisenberg picture, and obtain the sec- 
ond order differential equations 



U\ - 



—2y^ c — xy , 
-2xy [b^b + c^c) - 



ks 



(26) 



These equations imply that the first order derivatives of 
the operators 5^6, c^c and k^ change on a time scale of 
the order which is, for large N, much faster than 

the time scale on which the time evolution of b'^b, c^c 
and fca takes place. Since we are only interested in the 
time dependence of the phonon number b'^b, we can safely 
assume that the first order derivatives ^(6^6), ^(c'^c) 
and ^A:3 adapt adiabatically to the state of the system. 
Setting the right hand sides of the differential equations 
(|26|l equal to zero, we find 



c^c : 



b^ b and k-i 



2x 



b%. 



(27) 



y V 

Moreover, Eq. H27() implies, due to the conservation laws 
(113 and (123), that 



d_ 

dt ^ 



(&t6) = A(et,) ^ U. 



0. 



(28) 



This means that, after a short time of order the 
system reaches a quasi-stationary state with a constant 
ratio between the expectation value of the coherence k^ 
and the number of photons inside the cavity, respectively, 
compared to the total number of phonons in the system. 
To a very good approximation, these ratios remains con- 
stant throughout the whole coohng process. 



C. Cooling equations 

In the absence of spontaneous emission, at most a re- 
distribution of phonon energy can occur in the system 
0. Efficient cooling and the irreversible removal of en- 
ergy from the system requires dissipation. We therefore 
now consider the effect of dissipation in more detail. Us- 
ing the master equation JHJl we find 



dt 



dt' 



ix{S+b- S'b^). 

iy{S+b- S-b^) 
-^nks . 



K C , 



S-c^ 



(29) 

The leakage of photons through the cavity mirrors with 
decay rate k, not only decreases the number of photons 
in the cavity mode but also affects the size of ^3. 

In the previous subsection, we have seen that the sys- 
tem reaches a stationary state on a time scale of the 
order with c^c and fca being a multiple of Wb. Using 
Eq. we find 



[b'^b) 



dt 



(30) 




FIG. 2; Cooling of common vibrational modes obtained from 
a numerical solution of the master equation JSJ (solid line) in 
comparison to Eq. 13211 (dashed line) for g = 10~^ k, rjO, — 
5 • 10"" K, N = 10** and mo = 10^ (a). Cooling of individual 
phonon modes, obtained from a numerical solution of Eq. (21) 
in Ref. 0] under the assumption of negligible cavity photon 



numbers, for g 
mo = 10^ (b). 
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This shows that on the time scale we are interested in, 
namely a time scale of the order 1/Vn, both, the oper- 
ator c^c and the presence of a negative valued operator 
^3, provide effective cooling channels in the system. In- 
serting the results l|77jl and into Eq. ipn|l we obtain 
indeed 



Ttib'b) = 



■Kb^b. 



(31) 



If m denotes the phonon number expectation value {b^b)p 
and mo = rn{0), we finally obtain 



Too exp 



2 2 



Kt 



y^ 



(32) 



This describes exponential cooling of the atomic sample 



with a rate that can be as large as VNg and ^VNrjfl 
(cf. Eq. (UJ). The result ^ coincides with Eq. (17) 
in Ref. ,2] and is in good agreement with an exact nu- 
merical solution of the time evolution of the system (see 
Fig. [2Ia)). Spontaneous emission from excited atomic 
levels with decay rate F is too slow to contribute to the 
cooling of the system, since most of the particles remain 
in their ground states and k ^ F. 



D. Similarities between collective cooling and 
phase transitions 



Finally, we remark that the cooling process we describe 
here is not simply a redistribution of phonons into cavity 
photons and fcs, which then decay into the environment 
with the spontaneous decay rates k and ^k, respectively. 
Instead, the leakage of photons through the cavity mir- 
rors disturbs the equilibrium expressed by Eqs. (|23|l and 
(|25|l perturbing the otherwise conserved quantities 



and Q' = (fetfo) 



(33) 



The coupling to the environment causes a dynamical re- 
sponse in the system, which generates a transition from 
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a state with a fixed ratio of phonons to cavity photons 
(cf. Eq. (123) into a state with no phonons and photons 
in the setup. 

In atom-cavity systems, the time evolution of the sys- 
tem is, on a very short time scale, usually dominated 
by spontaneous emission. However, in the presence of 
many particles, the evolution is primarily governed by 
the Hamiltonian Hcomm, which drives the system into 
a quasi-stationary state characterised by the preserved 
quantities Q and Q' within a time of the order 1/7V. The 
presence of the spontaneous decay rate k, which scales 
as (cf. Eq. |p[J), disturbs this stationary state and 
causes the system to continuously assume new values for 
Q and Q' . Using Eqs. (EH) and (EH, we find 

While the relative size of the photon number (c^c) and 
the coherence (fcs) with respect to the number of phonons 
m remain constant, their sum, as accounted for by Q and 
Q' decrease exponentially in time due to the presence of 
dissipation, namely leakage of photons through the cavity 
mirrors. This is typical for a system undergoing a phase 
transition. 



V. COLLECTIVE COOLING OF INDIVIDUAL 
PHONON MODES 

There are many similarities but also many differences 
between the cooling of the common and the cooling of the 
individual phonon modes of particles. For example, one 
cannot simplify the Hamiltonian Hmd using the collec- 
tive lowering operator S", as we did in Eq. (|12|) . However, 
there is still a high symmetry in the system and all par- 
ticles are treated in exactly the same way. To take this 
into account we introduce the vector operators 

^^^^{""f)^ ^-^{^^) and c^{c} (35) 
with the usual scalar product, such that, for example, 

-h^-^yafh. (36) 

This notation allows us to write the Hamiltonian ((7)) as 

iJind = ^(a;5+ • 6 + yS:+ • c + H.c.) , (37) 

which is of a similar form as the Hamiltonian -ffcomm in 
Eq. 113). 

In analogy to Section llH Bl we proceed by introducing 
an effective annihilation operator for each particle i 
with 

ai='^{xb^ + yc) , a={a^] with [aj,a|] = l. (38) 

Then the Hamiltonian H37(l can be written as 

iJind = hzS+ -a + H.c. (39) 



Each particle couples individually to a new type of 
bosons, which are a superposition of a single phonon and 
a cavity photon. The number operator accounting for all 
bosons ai equals 

■ a= X: ^ ■ b + ^ ■ c + yx k:i \ , (40) 

where we defined the coherence 

fcg = 6t • c-f-6- . (41) 

While 6^ • h counts the total number of phonons in the 
system, the expectation value of • c = Nc^c. Again, 
the coherence describes a certain "symmetry" between 
phonons and photons in the setup and accounts for a 
continual conversion of the two types of bosons into each 
other. 



A. Conservation laws in the absence of dissipation 

Neglecting spontaneous emission, i.e. assuming k « 
and r « 0, and using the Hamiltonian H37(l . we find 

^^{h^ -h) = \x{S+ -b-Sr -6^), 

jih - iy{S+ -b-Sr ■b'<) + ixiS+ -c-Sr -cf). 

(42) 

In the derivation of these equations, we neglected the 
operators ajbj and (Ti&J with i ^ j since there are no 
interactions between particles and the phonons of other 
particles. From Eq. (|42|l we see that there are, as before, 
two conserved quantities in the system, namely 

Q ^ y^{b^ ■b) + ^{c^ ■c)-xy{k,), 
Q' ^ {b^-b) + ^{c^-c). (43) 

One is associated with the total number of bosonic parti- 
cles with annihilation operators ai, the other one counts 
the total number of phonons and photons in the setup 
and 

TtQ = TtQ'-0- (44) 

From the formal equivalence of the Hamiltonians 1)14(1 
and (|39|l pointed out in the previous subsection, one 
might have had expected that the expectation value 
{a} ■ a) is preserved in the time evolution of the system. 
However, this would only be the case if and com- 
mute with each other for i j and does not apply here. 

B. Quasi-stationary states 

To calculate the distribution of phonons, cavity pho- 
tons and coherence k^ in the system, which builds up 
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in the absence of spontaneous emission, we now consider 
the second derivatives of the operators b} ■ b, ■ c and 
fca in the Heisenberg picture. Neglecting again the popu- 
lation in the excited states, as accounted for by S_~^ ■ S_~ , 
we obtain 



and 



2x'^ X xy 



■(c^'c) — ~2y^ ■ c — xy , 



(45) 



Setting the right hand side of these equations equal to 
zero and taking the conservation laws (|44|l in the absence 
of dissipation into account, we find that the system pos- 
sesses a stationary state with 



y 



^(ctc) = 4.k 



dt 



dt' 



y 

3 = 0. 



(46) 



In contrast to the results in Section FlV Bl this equation 
describes a state with only a relatively small number of 
photons in the cavity mode compared to the total number 
of phonons in the setup. However, as we see in the next 
subsection, the presence of a negative coherence ^3 of 
the same order of magnitude as 6^ • 6 provides a cooling 
channel, which is sufficient to obtain a cooling rate of the 
same order of magnitude as the cavity decay rate k. 



niQ exp 



Kt 



(50) 



where to and toq now describe the total number of 
phonons in the system, as accounted for by 6^ • 6 at time 
t and t — Q, respectively. As in the case of the cooling 
of common vibrational modes, cooling rates of the same 
order of magnitude as the cavity decay rate k can be 
obtained, which can be as large as ^/N times the single 
particle coupling constants. 

However, there are some differences with respect to 
the case considered in Section Hvl In the present case, 
achieving such a high cooling rate and cooling the system 
to very low temperatures first requires the build up of a 
reasonably large coherence fcs, which does not exist in the 
absence of any laser driving. Suppose, there are no ini- 
tial correlations between the particles and their motional 
degrees of freedom. Then ^/ta = (cf. Eq. More- 
over, the second order derivative of fcs increases only on 
a time scale of order one (cf. Eq. (j^ ). Initially, the sys- 
tem is far away from the quasi-stationary state described 
in the previous subsection. However, once k^ reaches 
its equilibrium, the collective cooling process can begin. 
The coherence k^ needs no longer to be established; the 
system only has to adapt to the small changes of Q and 
Q' caused by the leakage of photons through the cavity 
mirrors. For a numerical solution of the time evolution 
of the system see Fig. |2Ib). 



C. Cooling equations 

As before, the effective removal of phonons from the 
setup requires leakage of photons through the cavity mir- 
rors. To take spontaneous emission into account we con- 
sider again the master equation JHJ and find 

^(6t .6) = ix{S+ -b-S- -b^), 

^(ct-c) = iNy{S+ -c-Sr Kc^ -c, 

^yfcg = iyiS+ -b-Sr ■ b^) + ix{S+ -c-ST ■ c') 
-^nks, (47) 

which implies 

Ttib'-b) = -^hct.c+^(ct.c)] 

As already mentioned above, the number of photons in 
the cavity remains relatively small and does not con- 
tribute considerably to the coohng process. Neglecting 
• c and using the results from the previous subsection, 
we obtain 

iih^-b) = -^^b^-b (49) 



VI. CONCLUSIONS 

The emergence of collective dynamics in a system of a 
large number of particles manifests itself in some macro- 
scopic features in the system behavior. As an example, 
we considered the problem of fast and efficient cooling 
of an assembly of N two-level particles trapped inside a 
leaky optical cavity. Results obtained here confirm those 
derived in Ref. 2]. The particles are excited by red- 
detuned laser fields. When the coupling constants are 
for all particles the same, a collective behavior emerges 
and the cooling rate can be as large as ^/N times the 
single-particle coupling constants. The generation of co- 
operative behavior of the N particles is crucial in the 
excitation step as well as in the de-excitation step. The 
collective states of the assembly of N particles then ex- 
perience a very strong coupling to the laser field as well 
as to the cavity mode. 

We have considered the case of particles sharing a com- 
mon phonon mode and the case of individual phonon 
modes for each particle. The two cases are similar and 
in both cases the collective cooling rate is very fast as 
compared to the single-particle cooling rate. The two 
cases are, however, different for the behavior of the co- 
herence k^: in the individual mode case, high cooling 
rates are achieved by first building up a reasonably large 
coherence k^, which does not exist in the absence of any 
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laser driving. Initially, the system is far away from the 
quasi-stationary state, where the phonon and the photon 
populations are balanced through the action. Only 
when ^3 reaches its equilibrium, the collective cooling 
process can begin and the system can adapt to the small 
changes of Q and Q' caused by the leakage of photons 
through the cavity mirrors. In the common mode case, 
there is no need for the "pre-cooling" phase for a build 
up of the coherence . 

The system considered in this paper may be a paradig- 
matic example for other applications of physical interest, 
such as systems undergoing continuous phase transitions, 
yet preserving some specific features during their evolu- 



tion. Systems presenting such a behavior might be of 
interest as well in biology. 
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